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Coherent optical control of polarization with a critical metasurface 
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We describe the mechanism by which a metamaterial surface can act as an ideal phase-controlled rotatable 
linear polarizer. With equal-power linearly polarized beams incident on each side of the surface, varying the 
relative phase rotates the polarization angles of the output beams, while keeping the polarization exactly linear. 
The explanation is based on coupled-mode theory and the idea of coherent perfect absorption into auxiliary 
polarization channels. The polarization-rotating behavior occurs at a critical point of the coupled-mode theory, 
which can be associated with the exceptional point of a parity-time (PT) symmetric effective Hamiltonian. 
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In photonics, optical loss is commonly regarded as an un¬ 
wanted nuisance. However, some recent advances have shown 
that loss is an interesting control parameter in its own right, 
and can be used to manipulate coherent light in useful ways. 
A case in point is the phenomenon of coherent perfect absorp¬ 
tion (CPA): when the loss in an optical structure is tuned to an 
appropriate (non-infinite) level, a specific incident wavefront 
is absorbed without scattering i Mllll . This is a generalization 
of the phenomenon of “critical coupling” fT|], and in a multi¬ 
channel system, like a metamaterial surface (“metasurface”) 
with waves incident from both sides, it provides a way to con¬ 
trol light with light without optical nonlinearity @-0]: varying 
part of the incident wavefront, such as the phase of one input 
beam, can switch the whole wavefront between perfect and 
near-zero absorption. Another example of optical loss as a 
control parameter comes from the field of PT symmetric op¬ 
tics, which deals with structures containing spatially-balanced 
regions of gain and loss lfl2l - (26ll . Such devices exhibit an 
unusual form of non-Hermitian symmetry-breaking El, and 
the “critical point” or “exceptional point” of PT symmetry¬ 
breaking fl4}l has been found to be associated with extraordi¬ 
nary behaviors like unidirectional invisibility lff5l - [l8ll23l - (26ll . 
Intriguingly, several links have been found between PT sym¬ 
metry and CPA. PT symmetric scatterers can simultaneously 
exhibit CPA and lasing fl9l[2lll22ll . and in some metasurfaces 
the occurrence of CPA can be mapped to the PT-breaking tran¬ 
sition of an effective Hamiltonian |27i|. 

The theory of CPA is agnostic about the nature of the 
loss I H, which can be some combination of Ohmic loss, 
fluorescence—or even radiation into other coherent channels 
lf28l l29tl . At first glance, treating radiative loss using the lan¬ 
guage of CPA may seem pointless, for the “absorption” of 
light from one input channel, and its complete transmission 
into another channel, occurs in so simple a system as a non¬ 
scattering waveguide. However, when there are multiple scat¬ 
tering channels, the CPA concept can provide an interesting 
method for the coherent manipulation of polarization |f3Cfl . 
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The control of polarization with pairs of coherent input 
beams has been explored in recent experiments by Mousavi et 
al. lITillm who showed that when two equal-power linearly 
polarized beams are incident on an appropriately-designed 
chiral metasurface (e.g. an arrray of asymmetrically-split wire 
rings), varying the relative phase tf> beween the beams can in¬ 
duce a complete rotation of each output beam’s polarization 
angle, with the output ellipticity varying by < 15°. Hence, 
the metasurface functions as a phase-controlled polarization 
rotator. To explain this, Mousavi et al. noted that, for a single 
input beam, the transmission is approximately circularly po¬ 
larized; say, left-circularly polarized (LCP). To explain the po¬ 
larization rotation, the reflection of the other input beam, inci¬ 
dent from the opposite side, must be right-circularly polarized 
(RCP). This implies, by mirror symmetry, that for each input 
beam the reflection and transmission have the same handed¬ 
ness. This seems counter-intuitive, for a chiral resonance with 
mirror symmetry along the propagation axis ought to emit to 
each side with opposite handedness. 

In this paper, we present a theoretical study of an ideal two- 
sided polarization-rotating metasurface, which reveals deep 
ties to the concepts of CPA and PT symmetry Il2l - 

[27h . The metasurface, differing from Refs. fill lulll, contains 
pairs of coupled resonators radiating into different linear po¬ 
larization channels. Using coupled-mode theory lf33l - [39ll . we 
show that when linearly polarized input beams are incident 
on each side, achieving perfect conversion to the other po¬ 
larization (i.e., 90° rotation) requires a specific balance be¬ 
tween radiative and non-radiative loss rates; this is analogous 
to ordinary CPA, which occurs at specific intrinsic loss levels 
lilll . For special “critical” choices of the frequency and loss 
parameters, the output beams become exactly linearly polar¬ 
ized for all values of the relative input phase <p , with polariza¬ 
tion angles varying with tf>. Under one-sided illumination, the 
reflected and transmitted beams have the same handedness, 
which results from interference between direct and indirect 
transmission processes. However, the critical metasurface is 
generally not a perfect circular polarizer under one-sided illu¬ 
mination, and the simple explanation involving phase-shifted 
RCP and LCP components f3l[ [32ll holds only in the limit 
where the radiative and inter-resonator near-field coupling 
rates dominate the dissipative loss rate, and the total loss is 
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negligible. Nonetheless, when two input beams are applied, 
the elliptically polarized reflection and transmission on each 
side combine to give zero total ellipticity. The critical points 
of the coupled mode theory, where this behavior occurs, are 
the PT-breaking points of a non-Hermitian effective Hamilto¬ 
nian, whose eigenvalues give the frequencies for the CPA-like 
perfect polarization conversion condition. (Similar mappings 
to PT symmetric Hamiltonians have previously been explored 
for CPA [[27], and for polarization conversion under one-sided 
transmission 0 .) Using other parameters choices, we can 
also switch the outputs between circular and linear polariza¬ 
tion by varying <p. We present full-wave simulation results 
verifying the predictions of the coupled-mode theory. 

Consider the plane metasurface depicted in Fig. [Ha), which 
is mirror-symmetric along the direction z normal to the plane. 
The metasurface is populated by pairs of coupled resonant 
modes, described by amplitudes q = [q x , q y ] T , which radiate 
in the x and y directions respectively. Plane waves are nor¬ 
mally incident on the metasurface from both directions, with 
wave amplitudes a = [a/ !+ , an -, a v +, « t -] T , where +(-) de¬ 
notes waves incident from the top (bottom) of the plane, and 
h(v) denotes the linear polarization component parallel to the 
x(y ) direction. Likewise, the waves leaving the metasurface 
are described by amplitudes b = [/?;,+, £>/,_, b v+ , b v -] T , with ± 
denoting waves exiting in the ±z directions. The coupling be¬ 
tween the metasurface resonances and the input/output waves 
is described by a set of coupled-mode equations Il33[ - [39ll : 


where 


-i£lq = K a 
K t q + Ca = b, 
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Here, 8^ (where fi e {x,y\) is the detuning of the operating 
frequency from the p-oricntcd mode, y* is the radiative scat¬ 
tering rate, y^ is the non-radiative dissipation rate, and k is 
the near-field coupling between the modes; all these parame¬ 
ters are real. The matrix K represents the radiative coupling 
between the metasurface and the input/output waves, while C 
represents the direct coupling between the waves. The forms 
of these matrices are constrained by the mirror symmetry of 
the metasurface, the definitions of energy and power in terms 
of the coupled-mode quantities, and optical reciprocity 0. 
The scattering matrix S, defined by S a — b, takes the form 

S = C + iK r n" I Ks|^; (3) 



FIG. 1: (Color online) (a) Schematic of the input/output wave ampli¬ 
tudes d and b, relative to the plane of the metasurface, (b)-(c) Proper¬ 
ties of the output beams emitted from the critical metasurface, versus 
relative phase 0 = arg(a^_/a* + ) of the equal-power input beams in¬ 
cident from each side. The metasurface is tuned to the k = T = 0.3, 
y d = 0.2 critical point, (b) Output power So - |b;,+ P + |b/,_| 2 , with 
the power of each input beam normalized to 1 , so that the total input 
power is 2. (c) Polarization angle t // = tan -1 ),^/^ i)/2. The outputs 
are linearly polarized (S 3 =0) for all 0. (e) Ellipticity parameter 
S 3 /So versus detuning <5 of a single input beam, under transmission 
from ±z to +z (solid curve), and reflection from ±z to ±z (dashes). 
At the critical point 8 = 0, neither reflected nor transmitted beam is 
exactly circularly polarized (S 3 /So 1 ). 


We now assume that the metasurface is designed so the res¬ 
onances have the same dissipation rates and frequencies: 

7 d x = 7 d y = y d , 8 X = 8 y = 8. (4) 

(The radiative scattering rates, however, can and will differ.) 
Furthermore, we consider purely /t-polarized incident illumi¬ 
nation, with input amplitudes <?/, = [«/,+, The outputs 

in each polarization are b/, = S A ah and b v — D «/,. Due to the 
mirror symmetry, S A has a symmetric eigenvector [1, l] r and 
an antisymmetric eigenvector [1,—l] r ; the latter, with eigen¬ 
value —1, corresponds to a node on the plane, and hence zero 
total loss. 

By varying the metasurface parameters and input ampli¬ 
tudes, it is possible to arrive at a situation where b /, = 0, 
i.e. all the /z-polarized incident light is re-emitted in the v po¬ 
larization and/or dissipated. This corresponds to CPA with 
v-polarized emission as one of the “absorption” channels (it 
would be “ordinary” CPA if y* = 0). For this to occur, the 
symmetric eigenvector of S A must have eigenvalue zero, and 
we can show that this occurs if and only if 


The sub-matrix S A (S v ) describes how light incident in the Ii(v) 
polarization scatters into the same polarization, while D de¬ 
scribes the cross-polarized scattering. 


y\ -yd = y s y + 7d = r, 

5 2 = zr-r 2 . 


(5a) 

(5b) 





















3 


We can satisfy Eq. (f5at by designing each resonator appropri¬ 
ately, as discussed below. Then if |/c| > T, perfect polariza¬ 
tion conversion can occur at two distinct frequencies. But if 
|/r| < T, Eq. (l5bl > cannot be satisfied for any 6. 

We can interpret the conditions (l5at— (l5bl) in terms of an 
effective Hamiltonian, via an argument from Ref. 10 . From 
Eq. ((3j. we write det(S A ) = det (H v - 5 ■ I) / [— det(fi)], where 


H v 


(i(jd-y' x ) k \ 

\ k Kyd + jy)) 


( 6 ) 


The eigenvalues of H v are the detunings for which det(S A ) = 
0. These detunings should be real, but H c is non-Hermitian. 
However, H t becomes PT symmetric [|7z] , with P = ct[ and T 
the complex conjugation operation, when Eq. ( l5al ) is satisfied. 
In that case, the eigenvalues of H v are ± V/r - T 2 , which are 
the solutions to Eq. ( l5bl >; then the eigenvalues are real for the 
PT-unbroken phase of H v , |/c| > T. In the PT-broken phase, 
|/r| < T, perfect polarization conversion cannot occur for any 
real 6. At k = +T, which are the critical points of Eq. ( l5bl > and 
the PT-breaking exceptional points jl4j] of H t , perfect polar¬ 
ization conversion occurs at only one detuning, 5 — 0. 

Now suppose the metasurface is tuned to one of the critical 
points, satisfying Eqs. (l5ab— (l5bl) with k — T, 6 - 0. For h- 
polarized inputs, the coupled-mode equations give 


bh± = 


“ (®h+ ' 


■ a h -) 


K± = -- 


1 /T 

2 


■yd 


r + yd 


(fl/ !+ + ah-) . 


(7) 
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The outputs have equal power, and third Stokes parameters 


Sf = — 21 m [b h± b*, ± ] 


r -yd \ah+\ 2 - \ah-\ 2 
r + yd 2 


(9) 


Hence, for equal-power inputs (|a ;,+| 2 = |a;,_| 2 ), both outputs 
are exactly linearly polarized (Sf - 0). 

Part of this result is easy to understand: the outputs are v- 
polarized for symmetric inputs (perfect polarization conver¬ 
sion), and /z-polarized for antisymmetric inputs (node on the 
plane). However, Eq. (0 goes further, and states that the out¬ 
put beams are linearly polarized for any choice of input beam 
phases. Varying the relative phase rotates the output beams’ 
polarization angles between [0, n/2]. This behavior is specific 
to the critical metasurface. (At the other critical point, k = - F , 
Eqs. 0 and 0 have opposite signs, and the phase shift rotates 
the polarization in the opposite direction. Thus, k 4 0 implies 
broken left-right symmetry on the metasurface.) 

The action of the critical metasurface as a rotatable linear 
polarizer is shown in Fig.[ljb)-(c). At cf) = arg(fl/,_/a/, + ) = 0, 
perfect polarization conversion occurs, and the power of each 
output beam reaches a minimum of (T - y,/)/(!' + y d ) due to 
dissipation. If the dissipation is weak, y d <sc I , the total loss 
is ~ 4y r i l \ and the polarization angle relative to the x axis is 
if/ as (n ±(f>)/2 for the outputs in the +2, direction. 

To better understand the critical behavior, we examine the 
reflection and transmission under one-sided /z-polarized illu¬ 
mination. As shown in Fig. [T]d), S 3 > 0 for reflection and 



FIG. 2: (Color online) (a),(b) Poincare sphere trajectories of the out¬ 
put amplitudes for the (a) +z and (b) —z output beams, for equal- 
power input beams with relative input phase <f> = arg(a;,_/«/,+). Ar¬ 
rows indicate the direction with increasing </>. The trajectories labeled 
I, II, III have k, 6 parameters given by the matching points in (c); the 
other parameters are T = 0.3 and y d = 0.05. (c) Phase diagram of 
the metasurface. The curve labeled L corresponds to the solutions 
of Eq. f5bl . for which /z-polarized input beams can be perfectly con¬ 
verted to the v polarization; these are also the real eigenvalues of 
the PT symmetric effective Hamiltonian ©. The curves labeled C+ 
and C correspond to solutions of Eq. OH- for which the inputs can 
be converted to circular polarization. The k < 0 part of the phase 
diagram is not shown here, but has a similar form. 


S 3 < 0 for transmission at the critical point, for either choice 
of input direction. From the definition of S 3 , this means that 
a single input beam produces reflected and transmitted beams 
with the same handedness. This is because the transmission 
is comprised of direct transmission of the linearly-polarized 
input, and re-emission from the metasurface resonators. At 
the critical point, the total transmission’s handedness is oppo¬ 
site to the re-emitted component, and the same as the reflected 
beam. Fig.fTJd) also shows that the transmission and reflection 
are not exactly circularly polarized for y t j 4- 0 , but these ellip- 
tically polarized components are nonetheless able to combine 
to form linearly polarized output beams. 

Away from the critical point, the metasurface ceases to act 
as a linear polarizer under equal-power incident beams, as 
each output beam becomes elliptically polarized. However, 
there is a remnant of the critical behavior, due to the winding 
behavior on the Poincare sphere. Right at the critical point, 
one cycle of <p winds each output amplitude along the equa¬ 
tor of the Poincare sphere, as depicted by the red loops in 
Fig. 0a)-(b). Away from the critical point, the trajectory no 
longer follows the equator exactly, but one complete cycle of 
</> still induces one winding of the longitudinal angle 2 if /—and 
hence a full rotation of the angle of the polarization ellipse’s 
semi-major axis, if/. This is true so long as the loops do not 
cross the poles (where the output beams become circularly po- 
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FIG. 3: (Color online) Design of the critical metasurface, (a) Re¬ 
flectance spectrum for a strip antenna of length a = 434.25 nm, 
width w = 50 nm, and thickness f = 30 nm, with incident beam 
linearly polarized parallel to the strip, (b) Reflectance spectrum for 
a split-ring resonator (SRR) of width w = 50 nm, main axis length 
b = 355 nm, arm length h = 105 nm, and thickness t = 30 nm, 
with incident light linearly polarized parallel to the main axis. Both 
antennas are free-standing; solid lines show full-wave simulation re¬ 
sults, while dashed lines show the least-squares fit to coupled-mode 
predictions, (c) Schematic of the metasurface unit cell, with the strip 
antenna aligned in x and the SRR aligned in y, separated by distance 
g. (d) Full-wave simulation results for S 3 /So, versus g and operat¬ 
ing frequency /, for the +2 output beam with symmetric equal-power 
input beams. The dashed line shows where the coupled-mode condi¬ 
tion Eq. ( 1 5h b is satisfied, using fitted coupled-mode parameters; the 
solid black line shows where S 3 = 0 in the simulation results. 

larized). It can be shown that the pole-crossings occur when 
6 2 = k 2 -F 2 ±2kJt 2 -y 2 . (10) 

For k > 0, the +(-) signs correspond to an RCP (LCP) +S out¬ 
put and an LCP (RCP) -z output. This results in the “phase 
diagram” shown in Fig. |2jc). The solutions to Eq. ( 1 Kil l lie 
along the curves labeled C±. In the region between these two 
curves, the longitudinal angle 2 i// undergoes a complete wind¬ 
ing with (f>. The behavior for k < 0 can be similarly deduced. 
When Eq. (ITOt is satisfied, varying </> switches the output be¬ 
tween linear polarization and circular polarization, as shown 
by the green and blue curves in Fig.[2ja)-(b). 

We have verified the above coupled-mode theory using nu¬ 
merical simulations of an exemplary plasmonic metasurface, 
shown schematically in Fig.[3jc). Each unit cell contains a sil¬ 
ver strip antenna radiating in the x direction, and a silver split¬ 
ring resonator (SRR) radiating in y. The cells are arranged in 
a square lattice with period d - 600 nm, and the entire meta¬ 
surface is free-standing in vacuum, which guarantees that no 
high-order diffraction mode exists below 500 THz. The di¬ 
electric function of silver is modeled by a Drude formula s m = 
£00 - f 2 p /(f 2 + iy P f), where f p = 2230 THz, y p = 5.09 THz, 
and Sc = 5. To extract the coupled-mode parameters for 
each resonator, we perform full-wave (finite-difference time- 
domain) numerical simulations of single-sided illumination 



FIG. 4: (Color online) Intensity So, ellipticity parameter S 3 /S 0 , and 
polarization angle 1 // for the +z output beam, versus the relative phase 
(f/ of the input beams. The metasurface is at the critical point indicated 
in Fig. [3jd). Solid curves show full-wave simulation results, and 
dashes show coupled-mode predictions using best-fit parameters. 


incident on each antenna separately, in the absence of the other 
antenna, with the appropriate linear polarization. The com¬ 
puted reflectance spectrum is fitted to the theoretical result 
R = y p 2 / [(/ - fp) 2 + (y p - y p ) 2 \ obtained from Eqs. (Ub>— 
(l2cl) in the k - 0 limit. Using the parameters stated in the cap¬ 
tion of Fig. [33 the resonators have equal resonance frequencies 
fx-fy- 206.3 THz and dissipation rates y d x = y^ = 1.3 THz, 
thus satisfying Eq. (0). Furthermore, the radiative decay rates 
are y s x = 19.8 THz and y x = 17.2 THz, satisfying Eq. (Bat . 

We then include both resonators in the metasurface, sep¬ 
arated by distance g, and perform another set of full-wave 
simulations. Varying g, with all other geometrical parame¬ 
ters fixed, alters the near-field coupling parameter k, as well 
as (weakly) the resonant frequency /o. Nonlinear fits of the 
simulation results to the coupled-mode theory give the func¬ 
tional relations k ~ 28.32 - 0.5lg + 0.0026g 2 , where k and g 
have units of THz and nm, and /o 200.48 + 0.073,g in THz. 
The fitted coupled-mode theory gives results for the various 
output Stokes parameters that agree well with the simulation 
results. For instance. Fig. [3jd) shows the simulation results of 
S 3 /So for one of the output beams, with a pair of symmetric 
input beams. The locus of S 3 =0 according to the simulation 
(black line) closely matches coupled-mode prediction (white 
dashes). Using Eq. d5bl> . we find that the critical point occurs 
at g w 22 nm, / ~ 202 THz. Fig.[4]shows the behavior at this 
critical point, which is in good agreement with the coupled¬ 
mode theory, particularly in the fact that S 3 « 0 for all </>. 

In summary, we have shown that the principles of CPA 
and PT symmetry can be used to design metasurfaces for 
coherently manipulating polarization using two input beams. 
This is an example of an emerging class of photonic devices 
that exploit the properties of “critical points” or “exceptional 
points” ll23l l26l l4lL l42ll . Here, the critical behavior corre¬ 
sponds to the output beams achieving exactly zero ellipticity. 
With other choices of metasurface parameters, other forms of 
polarization control can be achieved, such as switching be- 
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tween circular and linear polarization. Generalizations of the 
coupled-mode theory, such as having resonators which are not 
exclusively coupled to a single linear polarization channel, 
might be useful for designing other classes of polarization¬ 
controlling metasurfaces. 

M. Kang acknowledges supports from the National Nat¬ 
ural Science Foundation of China under Grant 11304226, 
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